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OQ . We use a new, conformally-invariant method of analysis to test incomplete null 

geodesies approaching the singularity in a model of an evaporating black hole for 
the possibility of extensions of the conformal metric. In general, a local conformal 
extension is possible from the future but not from the past. 
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1 Introduction 



In [7], Penrose developed his argument that, in the correct theory of space-time, initial 
singularities and final singularities would be different in character, with the Weyl tensor 
finite or even zero at initial singularities but no such constraint on final singularities. As 
a test case, he discussed the space-time singularities associated with Hawking's picture of 
an evaporating black hole (see also [8]). In this picture, there is a final singularity formed 
in the collapse which led to the black hole. At the endpoint of the evaporation of this 
black hole, this singularity vanishes but not before, at the last instant, being visible in 
the exterior and therefore being briefly an initial singularity. Penrose's argument then 
suggests that the singularity as seen from the past should be different in character, with 
respect to conformal properties, from the singularity as seen from the future. 

To test this suggestion, we need a convincing model of the metric for an evaporating 
black hole and we shall use the model of Hiscock [3], [I]. Essentially, this consists of 
matching Vaidya metrics to the Schwarzschild metric of the black hole and to flat space, 
with an ingoing flux of negative energy density which reduces the Schwarzschild mass to 
zero, and an outgoing flux of positive energy density which radiates the mass to infinity, 
leaving flat space to the future of an outgoing null hypersurface. Details will be given 
below, but note that the model is spherically-symmetric. 
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Once we have the model, we test its conformal properties using a method developed 
in [5]. Null geodesies are conformally-invariant as point sets, in the sense of being un- 
changed under conformal rescalings of the metric, and as a point set a null geodesic is also 
a null conformal geodesic (see e.g. [11] for the general definition of conformal geodesic). 
By regarding a null metric geodesic, say 7, as a null conformal geodesic, one can give 
a conformally-invariant definition of propagation of frames along 7 and can then give 
conformally-invariant definitions of boundedness of conformal curvature and its deriva- 
tives along 7, namely by requiring components of the (conformally-invariant) conformal 
curvature in these conformally-invariant frames to be bounded. In general, it is more 
natural to work with the calculus of tractors and tractor curvature (see [H]). For our 
purposes in this article this technology isn't needed, but we recall that the tractor curva- 
ture, which is the curvature of the tractor connection, has as nonzero components the Weyl 
tensor C abc d and the Cotton tensor, which is the derivative V[ a P fe i d of the Schouten tensor 
(defined below). Now suppose that 7 is incomplete as a metric null geodesic but that 
there is a conformal rescaling of the metric which allows an extension of the space-time 
so that 7 can be extended as a null geodesic of the rescaled metric. Then necessarily the 
conformal curvature (and the tractor curvature) will be bounded in the rescaled metric, so 
that the conformally-invariant conditions of boundedness will hold. (It is a recent result 
of one of us [5] that these necessary conditions are sufficient: given bounded derivatives 
up to order k + 1 of the tractor curvature, a C k local extension of the conformal structure 
will exist.) 

This method was used in [5] in a variety of particular cases: on radial null and space-like 
geodesies in the Schwarzschild metric, to find, as expected, that no conformal extension 
is possible through the singularities at r = or through space-like infinity io] on radial 
null geodesies of the conical spherical metric: 

ds 2 = dt 2 - dr 2 - a 2 r 2 (d6 2 + sin 2 6d(j) 2 ), a ^ 1, 

(so-called because it has a singularity at the origin due to a deficit in the solid angle, 
analogous to the deficit in the plane angle at the familiar conical singularity) to show 
that no conformal extension is possible through the singularity at r = 0; and along the 
matter flow-lines of the Lemaitre-Tolman-Bondi dust cosmologies to find conditions for 
conformal extendibility through the singularity. 

In this article, we use the method to analyse the singularity which is formed in the 
collapse and then evaporates with the black hole. This singularity has, in the Penrose 
diagram (see figure 1), an interior and a final point before it vanishes. Our findings may 
be summarised as follows: 

• along radial in or outgoing null geodesies meeting the singularity at an interior point, 
the physical affine parameter distance to the singularity is finite, and the Weyl spinor 
can never be made finite by rescaling; no conformal extension is possible (note that 
this includes regions of the Schwarzschild singularity, so in particular, as one would 
expect, no conformal extension is possible through there); 

• along radial ingoing null geodesies from the past to the final point of the singularity, 
the physical affine parameter distance is finite; the Weyl spinor vanishes by conti- 
nuity; the tractor curvature vanishes if one imposes Hiscock's condition ([3], [I]) of 
finite rates of particle creation; by imposing more conditions on the mass function 
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of the Vaidya interior metric, one may make more derivatives of the tractor curva- 
ture zero (and if all derivatives up to order k + 1 are bounded then there is a local 
C fc -extension of the conformal structure, j6]); if one imposes the decay rate on the 
mass function which follows naively from the Hawking mass-loss formula, then the 
tractor curvature is singular (and no local extension of the conformal structure is 
possible); 

• along radial outgoing null geodesies from the past to the final point of the singularity, 
the situation is similar to that with the conical spherical metric: there are choices of 
conformal factor which lead to a bounded Weyl spinor, but then the singularity is 
infinitely far off in the rescaled metric; all other choices lead to an unbounded Weyl 
tensor; either way no conformal extension is possible (it is also worth noting that 
the affine parameter distance to the singularity in the physical, unrescaled metric is 
infinite in this case); 

• along past-directed, radial ingoing null geodesies from the future to the final point 
of the singularity the physical affine parameter distance is finite, and of course the 
tractor curvature and its derivatives are zero; local conformal extension is possible. 

In summary, the conclusions accord with Penrose's hypothesis: a conformal extension is 
possible from the future but not from the past, with the exception that local extensions of 
low differentiability may be possible along ingoing null geodesies to the final point, given 
extra conditions on the rate at which the mass-function vanishes. A global extension from 
the past, through the whole singularity, is never possible. 

The plan of the article is as follows. We begin in the next section by giving more 
details of the evaporating black hole (EBH) models, introducing a Newman-Penrose null 
tetrad and calculating what we need from the spin-coefficient formalism (for which see e.g. 
j9] or [TO]). In section 3, we give the details of the method of testing for the possibility of 
a conformal extension, in terms of null conformal geodesies and the conformally-invariant 
propagation of spinors along them. In section 4, we collect the results of applying this 
test along radial in and outgoing null geodesies. 

2 The model 

As noted above, Hiscock's evaporating black hole (EBH) models [3], [I] are constructed 
by matching Vaidya metrics to the Schwarzschild metric of the black hole and to flat 
space. Start with a Schwarzschild metric, formed in collapse, and given to the future of a 
constant v surface S, say v = vg, in the usual coordinates (so v = t + r + 2m log(r — 2m) 
where m is the mass of the initial Schwarzschild metric). The details of the collapse are 
assumed to lie to the past of S. Pick a sphere S = {u = u±, v = v\) on S f = {v = v i > vg} 
and a spherically-symmetric time-like surface £ with past boundary at S. In the space- 
time region to the future of S' bounded by E, take the metric to be the ingoing Vaidya 
metric ([I]) below, with mass N(v), where N(v\) = m and N decreases to zero at vq (so 
N(vo) = 0) as negative mass flows into the hole. In the space-time region to the future of 
S" = {u = Ui} bounded by E, take the metric to be the outgoing Vaidya metric (j2J) below, 
with mass M(u), where M(ui) = m and M decreases to zero at u (so M(u ) = 0) as 
positive mass flows out to future-null-infinity X + . E has future boundary at the 2-sphere 
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Figure 1: Penrose diagram of the model EBH space-time; M, V + , V~ and Sch. are regions 
of Minkowski space, outgoing Vaidya, ingoing Vaidya and Schwarzschild respectively; V + 
meets V~ along the time-like hypersurface £, other matchings are along null hypersurfaces 
(V~ to Sch. along v = V\, V + to Sch. along u = U\, V + to M along u = uq); the 
singularity at r = (double line) becomes the origin of coordinates in M (dashed line) 
when the mass has dropped to zero (at v = Vq). 

S" — (u — u , v = v ). In the union of the regions to the past of S' and to the past of S", 
retain the Schwarzschild metric. 

There will necessarily be distributional curvature along £, and the choices of N(v), 
M(u), this distributional curvature and the history of £ are tied together by the Einstein 
equations. However, for our purposes, we don't need the details of this as long as there 
is a sphere (u = u ,v = Vq) on E as a future boundary at which N(vq) = = M{uq). 
Assume so, and then match to flat space across v = v and across u = uq. The Penrose 
diagram for this metric is given in figure 1. 

Note in particular that E must stop before the singularity is reached. This is because 
at the singularity r vanishes, while in the outgoing Vaidya metric the singularity at r = 
is in the past, corresponding to the past singularity of Schwarzschild, not the future one. 

Our interest is in the future singularity of the Schwarzschild metric at r = (the past 
one is in a part of Schwarzschild not relevant in collapse). Points of the future one with 
v < v i are just as in Schwarzschild, while points with v\ < v < Vq are as in Vaidya. Our 
main interest is in points with v = vq, and the approach to them along different radial 
null directions. When approached from the past like this, the last part of the approach is 
solely through the part of the metric which is ingoing Vaidya; approach from the future is 
through flat space; therefore we don't need to calculate with the part of the metric which 
is outgoing Vaidya, and so the details of E don't affect our calculation. 
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The ingoing Vaidya metric is 



ds 2 =[l- ) dv 2 - 2dvdr - r 2 (d6 2 + sin' 



r 

and, although we shan't need it, the outgoing Vaidya metric is 

ds 2 =(l- du 2 + 2dudr - r 2 {dd 2 + sin 2 9d0 2 ) (2) 

We choose the following Newman-Penrose (NP) tetrad for the ingoing Vaidya metric: 

D = d v + ~ (l - d r , A = -d r , (3) 

6 = -^=(d B + —„d d 



*V2 V sin^ 

so that £ a is outgoing into the future and n a is ingoing. Now calculate, using the NP 
formalism 

and 

where 

Note also that 



Dt = let , Dn a = -2en a , (4) 
A£ a = , An a = , (5) 

(6) 



N 



so that the surface r = 2N(v) consists of marginally outer-trapped surfaces of constant 
r, interpolating between the horizon of the initial Schwarzschild solution at r = 2m on 
v — v$ and the origin of flat space at v — Vq. This surface is a dynamical horizon in the 
language of pQ, but is diminishing in size as it has a flux of negative energy across it. 
The only nonzero curvature components are 

N' N 

$00 = -J" , ^2 = 3 ■ (8) 

so that 

$afe = ®00n a n b , ^ABCD = 6^/ 2 0(AObLc^D), (9) 

in terms of the spinor dyad underlying the chosen NP tetrad. 

Hiscock [3], [1] proposes imposing the condition N'(vq) = 0, as this is necessary and 
sufficient to give a finite rate of particle production at v — Vq. An alternative would be 
to take 

N(v) ~ (v - vf/\ (10) 

which is the dependence obtained if one assumes that the Hawking expression for the 
power radiated by an evaporating black hole holds up to the end. In this case, the Ricci 
tensor is singular at Vq. We'll return to this point below, but note here that, away from 
r = 0, Hiscock's condition has the effect of making the curvature continuous. 
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3 Conformal geodesies and the method of analysis 

A null conformal geodesic (see e.g. [IT] ) is a curve 7 with tangent vector v a and a covector 
b a defined along 7 satisfying the coupled system 

V v v a + 2{b c v c )v a = (11) 

V v b a -(b c v c )b a + ^(b c b c )v a = P ab v b (12) 

where P ab = & ab — Ag ab = — \R ab + j^Rgab is the Rho or Schouten tensor (here we 
use the conventions of [9] since we use the spin-coefficient formalism; other authors may 
have permutations of the signs here). By pip , a null conformal geodesic is a null metric 
geodesic with a particular scaling, dictated by a chosen solution of ffT2"j) . 
Define a propagation law along conformal geodesies by 

V v e a + {b c v c )e a + {b c e c )v a - {e c v c )b a = 0, (13) 

which as we shall see has a good transformation under conformal rescaling, then in par- 
ticular v a is carried along by this rule. For spinors the propagation ([TBI reduces to 

V v a A = -v AA 'b A , c a c (14) 

Under conformal rescaling 

g ab = tt 2 g ab (15) 
we claim that solutions of (TTT|) - (TT3|) transform according to 

v a = v a ,~b a = b a -T a ,e a = e a , 

where, as usual, T a = f2 -1 V a fi, and then (fT4"j) is invariant with a A = a A . Thus (1141) 
provides a conformally-invariant propagation of, say, a suitably normalised spinor dyad 
along 7, from which we construct a conformally-invariant test for boundedness of Weyl 
spinor components. 

For the normalisation, note that, given two solutions a± and a A of (|T4l) 

V v (e AB a A a2 ) = -{b c v°)e AB af af 

so that, if Q is a conformal factor with 

v a b a := v a (b a - T ) = (16) 

then 

- A R r\ AR 

is constant along the conformal geodesic. This rescaling of €ab corresponds to the rescaling 
(TT5]) of the metric and then (1111) becomes 

V v v a = 0, 

where V is the metric co variant derivative for g, so that v a is tangent to affinely-parametrised 
null metric geodesies for g. It is then easy to see that u a = Q 2 v a is affinely-parametrised 
for g. If s and s are affine parameters for u and v respectively then 

1 = v v s = ir 2 v u s = ir 2 ^. (17) 
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Note that s is a projective parameter along 7 in the usual terminology of conformal 
geodesies [TT], so that different choices of v a and b a subject to ffTTj) and ffT2]) lead to 
Mobius transformations in s. Equation (T3~6|) can be rewritten as 

and then the contraction of ffT2|) with w a gives 

V.V.IT 1 = n-\P ab v a v b ). (18) 

The method is now as follows: an affinely-parametrised null geodesic 7 of the space-time 
metric, with tangent u a , gives rise to a null conformal geodesic by solving (]12|) along 7 for 
6 a , with v a replaced by u a , and then adjusting the scaling by solving f[TTj) for t» a ; one can 
use ( 1171) to see if the null geodesic is incomplete for 5; if it is, one can use a spinor dyad 
solving (TH| to test the curvature components for boundedness on 7. If 7 is incomplete for 
s but with bounded curvature for g, then we expect to be able to extend the unphysical 
metric in a neighbourhood of 7 including a final segment [llj . [6]. This method therefore 
tests for extendibility of the conformal metric. 

We shall apply this method to ingoing and outgoing radial null geodesies in the ingoing 
Vaidya metric which encounter the singularity at r = 0, to test whether the conformal 
metric can ever be extended through this singularity. 

4 The analysis 

4.1 In-going null geodesies 

First we consider ingoing null geodesies, that is null geodesies parallel to n a . This geodesic 
vector field is affinely-parametrised for the physical metric (by so we take u a = n a . 
Note that the coordinate r is an affine parameter, so that the singularity is at a finite 
distance in the physical affine parameter, and the coordinate v is constant along such a 
7. We distinguish v = v 2 < v , an ingoing null geodesic meeting the singularity in its 
interior, from v = Vq, an ingoing null geodesic meeting the 'end' of the singularity. If 
v 2 < Vi, then we are testing the Schwarzschild metric for extendibility, while V\ < v 2 < v 
relates to the Vaidya metric. 

Expand b a in the NP tetrad as 

b a = X£ a + Yn a - Zm a - Zm a 

and expand (fT2l) using ([5]) and ([8]) to obtain the system 

AX = X 2 
AY = ZZ 
AZ = XZ. 

We won't need Y. Note that, from (116j) . 

n _1 Afi = b c u c = X 

from which we obtain Q, given X. 

This system of equations is homogeneous (since, by (jHJ), P a bn b = 0) so that the solutions 
are as in vacuum . For X and Z, use ([3]) to find two cases: 
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1. X = , Z = Zq for constant Zq, so = 1 w.l.o.g.; or 

2. X = (r + r ) _1 , Z = Z (r + Tq)" 1 for constants r and Zq, so f2 = (r + r ) _1 . 
Now for the spinor propagation ( 1141) . expand 



to find 



AC = 

At] = -Xr]-Z(. 
A basis of solutions in the two cases is given by 

1. A = o A + (770 + rZ ) L A ,I a = l a 

2. A = o A + (rj (r + r ) - Z )i A , I A = (r + r > A 

We wish to test the components of the Weyl spinor in these bases. From (Q in case 1 we 
find 

N 

ifj 2 = ^ 2 = - — 

which clearly diverges on the approach to the singularity at r = 0. In case 2, we have 

^2 = (r + r ) 2 ^ 2 

and, for any choice of tq, this still diverges if v = v 2 < vq, (which, by the remark above, 
includes the Schwarzschild singularity) but it will be zero on v = vq. To go further with 
the case v = vq, we calculate the other part of the tractor curvature, which is the Cotton 
tensor V[ a Pb] c and which will be conformally invariant at points where the Weyl tensor 
vanishes. For this Vaidya metric, it is given by 



V [Q P fe]c = V [a $ % 



with the aid of (jHJ) and (jHJ), which is 

N' 



(3£ [a n b] n c + n [a g b]c ) 



The component of this along I A I B O c O O is N'(v )/r 3 in case 1 and N'(v ){(r + 
ro) 2 /r 3 } in case 2, which diverges in both cases unless N'(vq) = 0. The vanishing of this, 
which is necessary for finiteness of the tractor curvature, is Hiscock's condition [3] , [1] for 
finite rates of particle creation at the end of the black-hole evaporation, so that there is a 
physical motivation for imposing it. If we do, then the whole of the tractor curvature is 
zero on this null geodesic (If we used (1101) instead, then the tractor curvature is singular 
and no extension is possible.) 

To test for the existence of a conformal extension, we need to consider higher deriva- 
tives of the tractor curvature. We next calculate the Bach tensor: 

B a b = 2V A ,V J?if> ABCD + ^ ABCD^ A'B' CD ■ 
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Recall that, under (|T5|) . B ab = Q 2 B ab . When v = vo, if N = N' = 0, the only nonzero 
term in B ab will be 

12 

0( A OBiCl<D)V a> Vi? iV(?;) 

2 .. 

= -^n a n b N , 

using ([3]). Now components of B ab in frames of either case diverge unless N" = 0, so 
there is no conformal extension possible unless N" = 0. Inductively, if N^ k '(vo) = for 
< k < n — 1 then the derivatives 

V ai • ■ • V afc C bcd 6 

will be zero for < k < n — 1, and the case k = n will be a conformally-invariant tensor, 
whose only nonzero term will be proportional to A^^ n V -3 ; so the (n — l)-th derivative of 
conformal (or tractor) curvature will be bounded if and only if N^ n \vo) = 0, and then a 
local C n_2 -extension of the conformal structure will exist, by [6]. 



4.2 Outgoing null geodesies 

Next we consider outgoing null geodesies, in the direction of l a . This vector field is not 
affinely-parametrised (see (j4])), so we introduce u a = Q 2 £ a which is. By (j4j), the condition 
on O is 

o = w u u a = e 2 D(e 2 £ a ), 

so that 

DO = -e6. (19) 
Again, we expand b a in the NP tetrad, to find for ffT2l the system 

DX + 2eX-ZZ = 

N' 

DY — 2eY — Y 2 = — (20) 
DZ-YZ = 0. (21) 

and for the spinor propagation 

D( + e( = -Y(-Zr) (22) 
Dr]-er] = 0. (23) 

For the conformal factor 

n~ l Dtt = b c t = Y, (24) 
Use (|24|) to solve (|2T|) and find Z = QZ Q , then (|20|) becomes a linear equation for Q^ 1 : 

N' 

D 2 n- X + 2e- 1 DODQ- 1 + = (25) 

which is equivalent to (jT8l) . Using ([7]), one solution is found to be Q = r -1 , corresponding 
to Y = p . The general solution can therefore be written as 



IT 1 = rF (26) 
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whereupon (|25p can be integrated to give 

DF = Ar~ 2 Q~ 2 (27) 

for constant A. 

From (1191) we may solve (1231) by 

7] = r/o©" 1 . 
Now use ([HID and (ED to convert ^ to 

A normalised spinor dyad is 

o A = eirV 
i A = wo A + e~h A 

where 

D(Q~ l QW) = -Z {1 2 9~ 2 . 
In this dyad, the nonzero Weyl curvature components are 

^ = tr 2 * 2 , (28) 

^ 3 = 3e- 1 Wn~ 1 ^2, (29) 
^ 4 = QW 2 Q- 2 ty 2 - (30) 

We need the behaviour of these along 7, where 

dr _ 1 N(v) 
dv 2 r ' 

which we need to solve for r(v). We distinguish two cases: 

• v — > v 2 < v : this includes the case of approach to the Schwarzschild singularity. 
Then 

r = {2N{v 2 )) 1 ' 2 {v 2 - v)V 2 (l + 0(v 2 - v)), 

so 

6 = (v 2 -v) l '\l + 0(v 2 -v)). 

From fr26|) - fr27|) . is therefore asymptotically 0(1) or 0((v 2 — v) 1 / 2 ) as v — > v 2 . 
From fl28|) and ([8]), we can't make ip 2 bounded. 

From the definition of G, the physical affine parameter s is obtained by solving 

£ = d, = e-, 

so that 

' d S = / 0" 2 ^;. (31) 



Now the distance to the singularity in the physical affine parameter is finite. 
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v ^ vq: now we need to make an assumption about the behaviour of N(v) subject 
to N(vq) = 0; we shall suppose that 

N(v) = n(v - v) k 

with k > 1, so that Hiscock's condition N'(vq) = holds, then along 7 we calculate 

r = 2n{v - v) k (l - 4nk{v Q - v) k ~ l + O((v - v) 2 ^)) . 
Solve (USD, using ®, to find 

9 = <* " -„)»-. ) (1 + °«"» " (32) 

From fl3Tl) and fl32l) . the integral for s diverges as v — > Vo, so that this singularity is 
infinitely far away along 7, as measured in the physical affine parameter. Note that 
this is different from the previous case of radial null geodesies ingoing towards this 
same singularity. 

From pi, (EHD and ©: 

F 2 N F 2 
r 2 



while, from ( 1271) 

dF A' 



ex P 1 Ti 7T? TJ-T (1 + 0((«b - «) )) (33) 



dv (vq — v) ik \4n(k — l)(v — v) k 1 
with A' = A(An 2 y 1 . 

We distinguish two cases. If A = then F is constant and ip2 is finite. By choosing 
Z = = VT, we ensure that the whole Weyl spinor is bounded (from (|29|) . (|30|) ). 
From (TTTj) however, choosing F = 1 we find that 

^1 = r - 2 0- 2 

and the right- hand- side of this is the same as in up to a nonzero factor. 

This equation has no solution which is bounded as v — > vq, so that the singularity 
is infinitely far off in the rescaled, unphysical metric - no conformal extension is 
possible. 

Now if A 7^ then we have to solve (133|) for F and again this has no solution which 
is bounded as v — > v , so that tp 2 cannot be made finite - in this case the distance 
to the singularity in the unphysical affine parameter is finite since 

^ = Q 2 Q~ 2 = r~ 2 F~ 2 Q- 2 = A^F- 2 ^, 
dv dv 

so that 

1 ° 2 

* = * + p, 
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which is bounded, but the Weyl curvature is singular so that no conformal extension 
is possible this way either. With Hiscock's condition, no local conformal extension 
is possible around the outgoing null geodesic to v — V\. 

It is straightforward to repeat the last calculation with (fTOj) instead. Now the 
physical affine parameter distance to the singularity is finite, but no choices lead to 
finite Weyl spinor, so that no local conformal extension is possible around this null 
geodesic. 
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